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Abstract. We review the status of mixing and lifetimes of fo-hadrons. We will show that AT/ AM, 
a s / and <p are better suited to search for new physics effects than AM alone, because of our poor 
knowledge of the decay constants. The theoretical precision in the determination of T^/M^ - which 
contains all information on AY/ AM, a s ( and <j> - can be tested directly by investigating the lifetimes 
of fo-hadrons, because both quantities rely on the same theoretical footing. In particular we will also 
present a numerical estimate for the lifetime of the S/,-baryon. 



INTRODUCTION - THEORETICAL TOOLS 



In this section we briefly discuss the principles of the calculation of physical quantities 
like lifetimes of b-hadrons x, the mass difference in the neutral B-meson system AM, 
the decay rate difference in the neutral B-meson system Ar, the semi-leptonic CP- 
asymmetry a s i and the mixing phase <j) . These quantities are currently measured at the 
B-factories and at the TeVatron (see e.g. [1]) and they will be measured at the LHC [2] 
or at a Super-B-factory [3] with high precision. They are defined as (see e.g. [4] for more 
details): 
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For the lifetimes one has to sum the decay rates into all possible final states X. The 
mixing stems from so-called box diagrams (see below). M\2 is the dispersive part 
(sensitive to heavy internal particles) and Ti2 is the absorptive part (sensitive to light 
internal particles) of these box diagrams. In the standard model T = 1 /T, A/12 an d Tn 
are given by the following diagrams (as an example we draw the diagrams for the B s - 
meson): 
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All these quantities are triggered by weak decays, in particular by the exchange of 
heavy W -bosons and the top-quark. Using the fact that these particles are much heavier 
than the b-quark (m t ,mw 3> m b ) one can integrate them out by performing an operator 
product expansion (OPE I), see e.g. [5] for a nice introduction. In the resulting effective 
theory the standard model diagrams are rewritten in a product of perturbative Wilson 
coefficients and new operators, they now look like that: 




c, u 



The vertices in the diagrams for T and T\2 are effective four-quark operators with 
AS = 1, while the vertex in the diagram for M\i is an effective four-quark operator with 
AB = 2. For M12 we have now already the final local operator, whose matrix element has 
to be determined with some non-perturbative QCD-method. 

As a next step we rewrite the expression for T in a form that is almost identical to the 
one of T\2. With the help of the optical theorem T can be rewritten (diagramatically: a 
mirror reflection on the right end of the decay diagram followed by all possible Wick 
contractions of the quark lines) in 
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The first term (=: To) corresponds to the decay of a free ^-quark, see e.g. [6, 7, 8] and 
references therein for some applications. This term gives the same contribution to all 
^-hadrons. The lifetime differences we are interested in will only appear in subleading 
terms of this expansion like the second diagram (=: r^), which looks very similar to 
the diagram for Y\2- Counting the mass dimensions of the external lines one can write 
formally an expansion of the total decay rate in inverse powers of the heavy quark mass 
m b : 

A A 2 A 3 

r = r +— r 1 + - I r2+- f r 3 + .... (4) 

m b m z b ml 

However the expressions for r ( - and Ti2 are still non-local, so we perform a second 
OPE (OPE II) using the fact that the ^-quark mass is heavier than the QCD scale 



(nib 3> Aqcd)- The OPE II is called the heavy quark expansion (HQE) [9, 10, 11, 12, 13, 
14, 15, 16]. The resulting diagrams for T 3 and T 12 look like the final diagram for M 12 : 



Now we are left with local four-quark operators (AB = for x and AB = 2 for T^). 
The non-perturbative matrix elements of these operators are expressed in terms of decay 
constants fs and bag parameters B. In the standard model one gets one operator for Myi, 
two independent operators for Tyi - including the operator that appears in Myi - and e.g. 
four operators for x{B + )/x{B ( i) 1 - in extensions of the standard model typically more 
operators arise. 

At this stage we would like to make some comments: 

• One can show that in the end no corrections of order I /nib survive in the total decay 
rate in Eq. (4). 

• r and T\2 are expected to be almost free from possible new physics contributions, 
since only light internal particles contribute, while M\2 might easily have large 
contributions from new physics effects. Since allowed new physics contributions to 
r and T\2 are smaller than the QCD uncertainties we neglect them in the following. 

• The OPE II seems to be theoretically less justified than the OPE I (mw/nib ~ 
17... 19 > 4... 10 ~ nib/AQCD), but the HQE can be tested directly by comparing 
experiment and theory for the lifetimes. 

• In all the diagrams shown in this section perturbative QCD-corrections have to 
be included! These corrections to the Wilson coefficients turned out to be quite 
sizeable. 

Summarizing one can state: the HQE represents a systematic expansion, which can in 
principle be tested by the lifetimes - in that sense it is not a model like the quark model. 



Besides testing our understanding of QCD and determining the standard model pa- 
rameters the search for new physics effects is a basic motivation for the study of the 
mixing quantities. Since M\2 is sensitive to heavy new internal particles one might start 
with AMj, which is proportional to /J B (see next section). Unfortunately fs s is hardly 
known. To visualise our current unsatisfactory knowledge of the precise value of the 
decay constant, we have taken some recent numerical values for fs s from the literature 
[17, 18, 19, 20, 21] and calculated the corresponding value of the mass difference AM S . 



This statements hold only at order l/m 3 b . 




INTRODUCTION - MOTIVATION 
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Depending on your favorite lattice collaboration you can arrive at theory predictions 
that are smaller, are equal or are higher than the experimental value for the mass 
difference AM,. This unfortunate situation might be called the decay constant problem. 
Here clearly more work has to be done to settle this ignorance and moreover the error 
estimates have to be done with much care. In the analysis in [4] we use the conservative 
estimate f Bs = 240 ± 40 MeV. 

In order to circumvent the decay constant problem one might try to determine the ratio 
AM s /AM d . Here the ratio \V t 2 jV t 2 d \ ■ fgBsJ '(/J^bJ arises. Although the ratio of the 
non-perturbative parameters is claimed to be theoretically better under control, one is 
still left with the uncertainty in the CKM elements, which is of the order of 40%. 
In the ratio r^/M^ the decay constant and the bag parameter B from M\2 cancel 
completely, schematically one gets 

We get a term that is completely free of any non-perturbative uncertainties (a) and a 
term that depends on the ratio of two bag parameters (b). If a > b and if the remaining 
uncertainties are under control then T\ilM\i might be an ideal quantity to search for new 
physics. Moreover the accuracy in the determination of this ratio can be tested directly 
via the lifetimes, which root on the same theoretical footing. 



STATE OF THE ART 



In this section we summarize the current status in the theoretical determination of the 
lifetimes of the &-hadrons and the mixing quantities. 



The mass difference - Myi 



Calculating the box diagram with internal top quarks one obtains 

M n , q = ^(V; q V tb ) 2 MiS (x t )B B jlM Bti f} B . (7) 

The Inami-Lim function So(x t = mf/M^) [22] is the result of the box diagram without 
any gluon corrections. The NLO QCD correction is parameterized by t)b ~ 0.84 [23]. 



The non-perturbative matrix element is parameterized by the bag parameter B and the 
decay constant f B . 



The decay rate difference - Tn 

The calculation of Ti2 is a little bit more involved since a second OPE has to be 
performed. Ti2 can be expanded as 



(8) 



The l/m£,-corrections (T^) were determined in [24] and they turned out to be quite 
sizeable. NLO QCD-corrections were done for the first time in [25], they also were quite 
large. At that time no lattice results were available for all appearing four-quark operators, 
so no real numerical prediction could be made. The first numerical estimate including 
NLO-QCD corrections and non-perturbative determinations of the appearing four-quark 
operators was given in [26] . Five years later the QCD-corrections were confirmed and 
also subleading CKM structures were included [27, 28]. Unfortunately it turned out that 
Ar is not well-behaved [29]. All corrections are unexpectedly large and they go in the 
same direction. This problem could be solved by introducing a new operator basis [4]. 
As an illustration of the improvement we show the expressions for T^/Mn in the old 
and the new basis: 
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Now the term that is completely free of any non-perturbative uncertainties is numerical 
dominant. Moreover the 1/m^-corrections became smaller and undesired cancellations 
are less pronounced. For more details we refer the reader to [4]. Currently also l/m/,- 
corrections for the subleading CKM structures in Ti2 [30] and 1 /m^-corrections for Ar, 
[31] are available - they are relatively small. 



Lifetimes 

The lifetime ratio of two Z?-hadrons can be written as 

i = i + ^r 2 + * (if + £r<" + ...) + ^ (rf + ...) + ... . (i i) 

T2 vanishes e.g. in x B +/xg d , ^B s /^B d and t s +/t e o but it survives in X\ b /xB d - The size- 
able NLO QCD-corrections to the lifetime ratios (T^) were determined in [32, 33]; 



1 /m£,-corrections (I4 ) and 1 /m^-corrections (ri ) were calculated in [34] - they are 
negligible for x B + / %B d and %bJ iB d , but they might be sizeable for Ta a / x Bd . 



NUMERICAL RESULTS 
Lifetimes 

The theoretically best investigated lifetime ratio is x B + / Xg d - One obtains [32, 33] 



r(B d ) 



= 1.063 ±0.027. (12) 



NLO-QCD corrections turned out to be important, while subleading l/m^-corrections 
are negligible. Some care has to be taken with the arising matrix elements of the 
four-quark operators: it turned out that the Wilson coefficients of the color-suppressed 
operators are numerically enhanced, see [32]. But the matrix elements of these operators 
are only knwon with large relative errors. Currently two determinations on the lattice are 
available [35, 36]. 

For TbJtbj large cancellations occur so the ratio is expected to be very close to one 
[24, 33] 

-7-4 = 1.00 ±0.01. (13) 
r(B d ) 

Predictions for the A b have to be taken with more care. In that case the NLO-QCD 
corrections are not complete and only preliminary lattice values [37] are available. A 
typical value quoted in the literature [38] is 

4^4 = 0.88 ±0.05. (14) 

The lifetime of the doubly heavy meson B c has been investigated e.g. in [39], but only 
in LO QCD. 

t(B c )lo = 0.5215;}| P s. 

In addition to the b-quark now also the c-quark can decay, giving rise to the biggest 
contribution to the total decay rate. 

An interesting quantity is the lifetime ratio of the S^-baryons, which was investigated in 
NLO-QCD in [32]. This quantity can in principle be determined as precise as x B + / x Bd 
(±3%). However, up to now the matrix elements for the E h baryons are not available. 
Assuming that the matrix elements for Z b are equal to the ones of A b we can give a rough 
estimate for the expected lifetime ratio. In order to get rid of unwanted s — > w-transitions 
we define (following [32]) 

i = r(Sfr) = r(z b ) - r(z b ^A b +x). (15) 



Using the preliminary lattice values [37] for the matrix elements of A b we obtain 



= 1-0.12±0.02±???, (16) 



where ??? stands for some unknown systematic errors. As a further approximation we 
equate t(S^) to t(A b ) - here similar cancellations arise as in TB s /TB d - , so we arrive at 
the following prediction 

= 0.88 ±0.02±???. (17) 



Mixing 

The mixing quantities have been investigated in detail in [4], numerically we obtain 



AM d = 0.53±0.18ps~ 1 , 


AM, = 19.3 ±6.7 ps" 1 , 
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Ar d /AM d = (52.6l|i; 5 8 )-10- 4 , 


AT S /AM S = (49.7 ±9.4) • 10~ 4 , 
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to = -0.09l!8;jgf, 


5 = (4.2±1.4)-1O- 3 , 
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a d fs = (-4.8 + _[° 2 )-10-\ 


a s fs = (2.06 ±0.57)- 10~ 5 . 


(23) 



The predictions for Ar^ and AT d /T d are obtained [4] under the assumption that there 
are no new physics contributions in AM d . From this list one sees the strong suppression 
of and a s i in the standard model. 



EXPERIMENTAL STATUS 
Lifetimes 

The Heavy Flavor Averaging Group quotes [40] the following numbers 

1.071 ±0.009, 4^4 = 0.939 ± 0.021, (24) 
?{B d ) 

0.921±0.036, t(B c ) =0.463 ±0.071 ps. (25) 

From the ratio x B + / %B d it can be seen that the HQE works very well. TbJ %B d is about 
2.9 o below 1, here more precise numbers are needed, to see whether there might 
be some interesting effects. The situation for the A^-baryon is not settled yet. First 
several theoretical improvements have to be included, second there are two different 
experimental numbers on the market [41, 42]. For B c the number lies in the right 



T(Afc) 

<B d ) 



ball park, but here also a full NLO-QCD calculation would be desireable to make the 
comparison more quantitaive. Finally we are waiting for a first result for the lifetimes of 
the S^-baryons. 



Mixing 

The mass differences have been measured with great precision at LEP, TeVatron and 
the B factories [43, 44, 45, 40] 

AM d = 0.507 iO.OOSps- 1 , (26) 
AM S = 17.77 i0.10i0.07ps" 1 . (27) 

Due to the uncertainties in the decay constants, theory will not be able to achieve a 
similiar accuracy in the foreseeable future. 

For the remaining mixing quantities a sh AT and <j> we do not have measurements yet, 

but very interesting bounds: 

In [46, 47] the dimuon- asymmetry was determined 

a d = 0.582flJ + 0.418a^, (28) 

a D0 = (-5.3i2.5il.8)-10~ 3 , (29) 

a^ DF = (i8.0i9.0i6.8)-10~ 3 . (30) 
The semileptonic CP asymmetry was also measured directly in [48]. 

a s sl = (1.23i0.97i0.17)-10~ 2 . (31) 

Here more precise numbers are needed, because a clear deviation from the small stan- 
dard model value would be be an unambiguous sign for new physics ! The same argument 
holds for the phase <j> s , while a clean measurement of Ar is probably best exploited by 
comparing experiment and theory for AT /AM. 

AT and <j> s have been determined from an angular analysis in the decay B s — > J/ij/Q: In 
the untagged analysis from DO [49] the following values were obtained 

0, = 0.79i0.56±g : J?, (32) 
Ar = 0.17 i 0.09 i 0.02 ps" 1 . (33) 

One has to keep in mind the 4-fold ambiguity in <j> s : with (j) s also — <j) s and 7t±(j) s are 
solutions ! CDF obtained from the untagged analysis [50] 

Ar = 0.076+g;Jgi 0.006 ps -1 . (34) 

and no bound on (j) s . 

CDF also performed a tagged analysis [51] and obtains confidence regions in the <j> s — 
Ar-plane, which differ about 1.5 cr from the SM prediction. If they fix |Ti2| to the SM 
value obtained in [4] they get 

-<j) s G [0.24, 1.36] U [1.78,2.90]. (35) 
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FIGURE 1. Experimental bounds in the complex A,-plane (state: end of 2006). The bound from AM S 
is given by the red (dark-grey) ring around the origin. The bound from AT S /AM S is given by the yellow 
(light-grey) region and the bound from a S f S is given by the light-blue (grey) region. The angle S A can be 
extracted from AF S (solid lines) with a four fold ambiguity - one bound coincides with the x-axis! - or 
from the angular analysis in B s — > J/y(j) (dashed line). If the standard model is valid all bounds should 
coincide in the point (1,0). The current experimental situation shows a small deviation, which might 
become significant, if the experimental uncertainties in AF S , a", and <j> s will go down in near future. 



NEW PHYSICS MODELS 

In the literature many new physics models are applied to the mixing sectors, e.g. [52, 53] 
and references in [4] . In [4] we have presented a model independent way to determine 
new physics effects in the mixing sector. We assume that new physics does not alter F\2 
- at least not more than the intrinsic QCD uncertainties, but it might have a considerable 
effect on Mi 2. Therefore we write 

r 12 = rf 2 M M 12 =M? 2 M -A (36) 

By comparing experiment and theory for the different mixing observables we get bounds 
in the complex A-plane, see [4] . Taking the solution for (f) s from the untagged DO analysis 
in the 4th quadrant - which corresponds to a certain choice of the strong phases in the 
decay B s — > // - and the data that were available at the end of 2006, we obtained in 
[4] a 2 a deviation from the standard model, see Fig. (1). A new analysis is currently in 
progress. 

Note added: (taken from [54]) There is sometimes a confusion between the mixing 
phases fi s and Q s , which we would like to adress here. Both numbers are expected 
to be small in the standard model - <j) s = (0.24 ±0.04)° and 2/3 5 = (2.2±0.6)°(= 



(0.04±0.01)rad), but in view of the high future experimental precisions - in particular 
at LHCb [2]- a clear distinction might be useful. 

2ft := -arg[(V tb V*) 2 /(V cb V* s ) 2 ] is the phase which appears in b — > ccs decays of neutral 
B-mesons taking possible mixing into account, so e.g. in the case B s — > J/\\f+ty. 
(V tb V* s ) comes from the mixing (due to M12) and (V c bV* s ) comes from the ratio 
of b — > ccs decay and b — > ccs amplitudes. Sometimes ft is approximated as 2ft ps 
— arg[(y^y f *) 2 ] ps — arg[(V ; *) 2 ] - the error due to this approximation is on the per mille 
level. 

<j> s := argfM^/r^] is the phase that appears e.g. in a^ s . In Myi we have again (V tb V t * s ) 2 , 

while we have a linear combination of (V cb V* s ) 2 , V cb V* s V ub V* s and (V^V^* ) 2 in Ti 2 . 
Neglecting the latter two contributions - which is not justified - would yield the phase 
2ft. 

New physics alters the phase —2ft to A — 2ft and the phase <j) s to A + (j) s . If the new 
physics contribution is sizeable, then in both cases only A survives, since the standard 
model phases are very small. 

In the tagged analysis CDF [51] introduces the phase 2ft for which the following 
relation to the notion in [4] holds -2ft := A - 2ft 5M . 

OUTLOOK 

In this talk we have summarized the current theoretical status of the lifetimes of b- 
hadrons and the mixing quantities. Our main strategy for finding new physics in these 
quantities is the following: New physics is expected to have the biggest effects in Mi 2, 
but due to the decay constant problem the quantity that comes first in mind - AM - seems 
to be not the best choice. We have argued that r 12 /M 12 is theoretically very well under 
control. Therefore our first choice are the quantities AT /AM, a s i and 0. Moreover the 
theoretical precision in the determination of T\2 can be tested directly by investigating 
the lifetimes of ^-hadrons, because both quantities rely on the same theoretical footing. 
We conclude with a subjective wish-list for theory and experiment: 

• Perturbative calulations: 

- NLO-QCD corrections for t(B c ) 

- complete NLO-QCD corrections for z(A b ) 

- r^forr^ 

- rfhorTn 

• Non perturbative calculations: 

- matrix elements for t(5 + ) / x(B d ) 

- matrix elements for x{A b ) and t(S^) 

- precise and relieable values for the decay constants 

- 1 /m-operators for Tn, a first step in that direction has been performed in [55] 

• Experiment: (ranked) 

1) Precise values for a s i and <j> s 

2) Precise values for t(B s ) / t(5j) and t(A&) and Ar — a first value for t(E&) 
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